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ABSTRACT 

It is shown that the acceleration of the 
solar wind particles at the solar wind 
terminal shock is capable of providing 
the total flux, spectrum and radial gra- 
dients of the low-energy protons close 
to ones observed in the interplanetary 
space. 


1 . Introduction. The solar wind deceleration by the in- 
terstellar medium may result in the existence of the solar 
wind terminal shock [l]. In this case a certain fraction of 
thermal particles after being heated at the shock would ob- 
tain enough energy to be injected to the regular acceleration 
process [2-5] . An analytical solution for the spectrum in the 
frame ef a simplified model that includes particle accelera- 
tion at the shock front and adiabatic cooling inside the 
stellar wind cavity has been derived in [6j. 

A more realistic model is discussed in the present pa- 
per in order to estimate whether this process might occur 
the main source of the low-energy cosmic rays in the helio- 
sphere . 

2. Model. Let the solar wind velocity u(r) in a steady 
state spherically symmetric model be equal to u^sconst at 
r < R and u 2 (fi/r) 2 at r > R, where R is the solar wind ter- 
minal shock radius, & = u^/u 2 - the compression ratio of 
the shock. Cosmic ray propagation in the diffusive approxi- 
mation is governed by the transport equations 
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where n(r,p) is the particle number density per unit inter- 
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val of momentum p, ^ (r,p) - cosmic ray diffusion coeffi- 
cient. At r < R and r > R (subscripts 1 and 2, respectively) 
it is given by 


The matc hing conditions are that the particles current 
and the number density should be continious at the shock. 

The other conditions are the absence of source or sink at 
r=0 (r-n(r,p)-> 0 as r -*• 0) and a given spectrum at infinity 
(n(r,p>-* N^p) as r-*©° ). 

Integrating (1) with respect to r from 0 to R-0 and 
using the analytical solution to (1) at r > R, the matching 
condition might be reduced to the following: 
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where f(p) ~ L~ ^ /( ^A) 

Hence, the problem is to solve the equation (1) at 
r ^ R with respect to the boundary condition (2). The equa- 
tion (1) is valid at particle momentum p higher than some 
limiting value, say p Q . One more condition is necessary to 
make the system (1),(2) being completed, namely n(R,p 0 )=n 0 . 
This condition relates the injection rate to the spectrum 
amplitude p=p p . 

It has been shown [7-9] , that the shock acceleration 
is efficient enough to transfer a significant fraction 
( ~ 50 %) of the inflowing plasma kinetic energy to the acce- 
lerated particles. The accelerated particles spectrum takes 
on the power law form with the index equal to that one given 
by the linear theory at particles velocity v or 5 ^ [9]* 

The corresponding values for the solar wind particles are 
T ^ 10 keV or p o ~ 5*10 ni c. Particles of energy correspon- 
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ding to the cutoff velocity diffuse far upstream and suffer 
catastrophic adiabatic cooling. 

3. Results and Discussion. The results of the numerical 
calculations are presented in Figure in comparison with the 
observational datafio]. The curve labelled Y' is the differen- 
tial proton flux near the shock, "2"- that one at the orbit of 

the Barth. The following parameters were used* ^=5*10^ cm/s, 

•14 , 1 

B=7,5»10 cm, o =4. The solar wind energy density near the 

shook E k has bean taken equal to 5 eV/cm 3 . The value of n„ 
has been defined so, that accelerated particles (T ^10 keV) 
contain 1C$ of the kinetic energy density E^. The diffusion 
coefficients are K =6,5«10 22 (p/mpC) 0, 2 (r/R) 0 * 2 cm 2 .s”^, 
ie z <<^ ^e i . The diffusion 
coefficient in the out- 
side vicinity of the 
shock might be signifi- 
cantly reduced due to the 
magnetic field compres- 
sion and convection of 
the turbulence excited 
by the accelerated par- 
ticles. This leads to the 
strong modulation of the 
galactic proton flux at 
T ^ 10 MeV. Radial gra- 
dient n r takes 

on a value of 15 %/AU at 
T— 1 0 keV and 5%/AU at 
T * (l-IO)MeV. 

The comparison of 
the results with the ob- 
servations shows that 
under the above assump- 
tions the acceleration 
rate (source (£< ) is ca- 
pable of providing the 
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registered intensity of the low-energy particles in the in- 
terplanetary space. In the energy range T 100 KeV the solar 
cosmic rays might contribute significantly to the total in- 
tensity, especially during the maximum of solar activity [ll]. 

This would cause the complicated variations of the low-ener- 
gy cosmic ray intensity: in the energy range T < 100 KeV the 
intensity is correlated with the solar activity level [10] j 
at the energy T=6 MeV/nucl. the anticorrelation takes place 
[ll]. The relative contribution of the source increases 
with the heliocentric distance r. 

The cutoff energy for the accelerated particles spectrum 
is determined by X< (T mfly ) - Ru^ ; on the other hand, ^ 

10 MeV/nucl. as is observed for different nuclear species. 

That shpuld be consistent with our model if Je. is a func- 
tion of T per nucleon only. The verification of this sugges- 
tion should be based on the treatment of the cosmic ray che- 
mical composition. 
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